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We apply a newly proposed Monte Carlo method, the Wang-Landau algorithm, to the study 
of the three-dimensional antiferromagnetic g-state Potts models on a simple cubic lattice. We 
systematically study the phase transition of the models with 5=3, 4, 5 and 6. We obtain the finite- 
temperature phase transition for q= 3 and 4, whereas the transition temperature is down to zero 
for q=5. For g=6 there exists no order for all the temperatures. We also study the ground-state 
properties. The size-dependence of the ground-state entropy is investigated. We find that the 
ground-state entropy is larger than the contribution from the typical configurations of the broken- 
sublattice-symmetry state for q = 3. The same situations are found for q — 4, 5 and 6. 

PACS numbers: 75.10.Hk, 05.50.+q, 05.10.Ln 



I. INTRODUCTION 



The g-state Potts model has various interesting prop- 
erties to study ■ The order of the phase transition 
of the Potts model depends on the spatial dimensional- 
ity and the number of states, q. The phase transitions 
of the antiferromagnetic (AF) Potts models are more 
complex than those of the ferromagnetic Potts models. 
Here, we focus on the three-dimensional (3D) AF g-state 
Potts models. Banavar et al. Q studied the AF 3- and 
4-state Potts models by use of the Monte Carlo sim- 
ulation. Ono J^j pointed out an appropriate choice of 
the vector order parameter for the AF Potts models on 
a bipartite lattice. The phase transition of the AF 3- 
state Potts model is considered to belong to the XY 
universality class HQ. The low-temperature phase of 
the AF 3-state Potts model was shown to be the broken- 
sublattice-symmetry (BSS) state, but it is still an open 
question whether there exists the rotationally symmetric 
state or not below the 2nd order phase transition point 
of the XY universality class Compared to the 

AF 3-state Potts model, there have been not so many 
systematic studies on the AF 4-state and higher-state 
Potts models in three dimensions. Recently, the AF 4- 
state Potts model was studied by Itakura fl2|l ; there is a 
finite-temperature phase transition, and the critical phe- 
nomena may belong to the Heisenberg universality class 
if the phase transition is of 2nd order. However, a pos- 
sibility of the 1st order transition was also argued [[L2[ . 
Another source of interest in AF g-state Potts models 
is that it has the nonzero ground-state entropy without 
frustration. Nonzero ground-state entropy, So 7^ 0, is an 
important subject in statistical mechanics. One physi- 
cal example is provided by ice. The AF Potts model is 
served as a useful model for the study of the nonzero 
ground-state entropy, and the nonzero ground-state en- 
tropy of the two-dimensional (2D) AF Potts model was 
extensively studied by Shrock and Tsai Jk| . 

The Monte Carlo simulation is regarded as a standard 



tool for studying statistical mechanical properties p4[ . 
The Monte Carlo study using a conventional Metropo- 
lis algorithm sometimes suffers from the problem of slow 
dynamics, or the long time scale problem. In the study 
of the phase transition, the long time scale due to the 
critical slowing down causes the problem of slow equili- 
bration. For the simulational study of the AF Potts mod- 
els, the cluster algorithm |1| has been used §^,|,|^,|2) 
to overcome the critical slowing down. As an extension 
of the cluster algorithm, Tomita and Okabe |l6| recently 
proposed a new cluster algorithm of tuning the critical 
point automatically, and applied it to the study of 2D 
Potts models. The extended ensemble method is another 
type of attempts to overcome the slow dynamics. The 
multicanonical method |l7| , |l8|| , the simulated tempering 
]l9t , the exchange Monte Carlo method j^fj] (or the mul- 
tiple Markov chain method [0), the broad histogram 
method |22j , and the flat histogram method |^3[ are ex- 
amples of the extended ensemble method. Quite recently, 
Wang and Landau Q proposed an efficient algorithm to 
calculate the density of states with high accuracy. In the 
present paper we use the Wang-Landau algorithm |24|] 
to study the 3D AF Potts models. We study not only 
the phase transition of the AF Potts models but also the 
ground-state properties. The Wang-Landau algorithm is 
quite effective for this purpose because we can calculate 
the density of states with high accuracy. We make a ran- 
dom walk in the whole energy space, and the algorithm 
is appropriate for calculating the ground-state entropy. 



We organize the rest of the paper as follows. In Sec. 

II, we describe the model and the vector order parame- 
ter for the AF Potts model. We also briefly explain the 
simulation method, the Wang-Landau algorithm. In Sec. 

III, we study the phase transition of the AF 3-, 4-, 5- 
and 6-state Potts models on a simple cubic lattice. The 
ground-state properties of the 3-, 4-, 5- and 6-state Potts 
models are studied in Sec. IV. The summary and discus- 
sions are given in Sec. V. 
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II. MODEL AND SIMULATION METHOD 

We deal with the AF Potts model defined by the 
Hamiltonian 



(1) 



<»j> 



where the Potts variable di takes the value 1,2, ■■■,q, 
and the summation is taken over the nearest-neighbor 
pairs of sites on a simple cubic lattice. 

The ferromagnetic 3-state Potts model can be trans- 
formed into the Z3 clock model, and the order parame- 
ter is well represented by a two-dimensional vector. For 
the AF Potts model on a bipartite lattice, the sublattice 
structure should be taken into account [§,^5); that is, 
the staggered magnetization will be treated. Using three 
components of the staggered magnetization 



E 



E 

jeB 



N, 



(2) 



where A and B denote two sublattices, we define the vec- 
tor order parameter as follows: 



M 



mi ex + m 2 e 2 + m 3 e 3 . 



(3) 



Here, denotes the unit vectors in two dimensions di- 
recting 120 degrees apart from each other, and 



ex + e 2 + e 3 = 0. 



(4) 



Then, the square of the vector order parameter simply 
becomes 



M = - (mi + m\ + ml). 



(5) 



We can extend the above argument to higher-state Potts 
models ||. In general the order parameter of g-state 
Potts model is well described by the (q — 1) dimensional 
vector order parameter. We consider the unit vectors in 
(<7 — 1) dimensions pointing to q directions as in the case 
of the 3-state Potts model; the sum of the unit vectors 
is set to be zero as in Eq. (Q). Using each component of 
the staggered magnetization, Eq. (0), we can obtain the 
generalized expression for Eq. (ph , 



M 



9 



9-1 



(6) 



We should note that for the 3-state Potts model the vec- 
tor order parameter space spans a hexagonal region in a 
two-dimensional space, and the maximum value of |M| 
is v3/2. For the 4-state Potts model the order param- 
eter takes the values within the area shown i n Fig . 3 of 
Ref. |l2), and the maximum value of |M| is y/2/3. 

We briefly describe the Wang-Landau algorithm 
This algorithm is similar to Lee's version of the mul- 
ticanonical method (entropic sampling) [fill, the broad 



histogram method p2[ , and the flat histogram method 
p3f ; but the Wang-Landau algorithm has advantage that 
it can estimate the density of states efficiently even for 
large systems. The idea of the Wang-Landau algorithm 
p4j is that we make a random walk in energy space based 
on the transition probability from energy level E\ to E 2 ; 



p{E\ — * E2) = min 



gCgi) 
g{E 2 y 



(7) 



where g(E) is the density of states. Since the exact form 
of g{E) is not known a priori, we determine g(E) itera- 
tively; g(E) is modified by 



In g(E) ^ In g(E)+ In fr, 



(8) 



every time the state is visited. The modification factor 
fi is gradually reduced to 1 by checking the "flatness" of 
the energy histogram; the histogram for all possible E is 
not less than some value of the average histogram, say, 
80%. 

We simulate the AF q-state Potts model (<?=3, 4, 5, 6) 
on a simple cubic lattice by using the Wang-Landau algo- 
rithm Q . We impose the periodic boundary conditions 
and the linear sizes are L =8, 10, 12, 14 and 16. For the 
modification factor f i: we start with / = e k with k=l 
or some positive integer, and /j+x = \fTi'i the final value 
of fi is chosen as 10 _s , which is the same as Ref. [ p4[ . 
We calculate the density of states g(E), and measure the 
physical quantities of interest as a function of E. Then, 
the canonical average of the physical quantity Q at the 
inverse temperature (3 = l/fc^T is calculated thorough 
the standard relation 



(Qh = 



jQ(E)g(E)e- 



P E dE 



Jg(E)e-f 3E dE 



(9) 



In the actual calculation, the relative density of states, 
g(Ei)/g(E 2 ), is directly obtained. In terms of the en- 
tropy (in units of fcs), S(E) — lng(E), the entropy dif- 
ference, S(Ei) — S(E 2 ), is directly measured. Imposing 
the constraint 



(10) 



we can determine the absolute value of g(E). Here, 
N(= L 3 ) is the number of lattice sites, that is, the num- 
ber of Potts spins. For the AF Potts model on a simple 
cubic lattice, Eq. (|l|), the energy E takes the value from 
to 37V in units of J. The state with the highest en- 
ergy E = 3N is nothing but the ferromagnetic ground 
state, and the degeneracy of the ferromagnetic ground 
states is q. Therefore, we may check the accuracy of the 
calculation by confirming <?(3A) = q. 
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III. PHASE TRANSITIONS OF AF POTTS 
MODEL 

First, we study the phase transition of the AF Potts 
model with <?=3, 4, 5 and 6. Let us start with showing 
the data for q — 3 in order to make a comparison with 
higher-state models, although the phase transition of the 
3-state Potts model has been studied extensively [p|-^0|. 
The temperature dependence of (M 2 ) for q=3 is shown 
in Fig. 0(a). The data for L=8, 10, 12, 14 and 16 are 
plotted by the dot-dashed, dashed, short-dashed, dotted 
and solid lines, respectively. The temperature is repre- 
sented in units of J/ks- All the measurements are done 
for 4 independent runs, and the average is taken over 4 
samples. The normalized fourth-order cumulant of the 
magnetization, the Binder parameter 



1 



1 - 



q-1 (M 4 ) 



9 - 



)■ 



(11) 



is plotted in Fig. |l|(b). The normalization factors in 
Eq. (|TT|) are chosen such that 



1 for T = 
for T = oo 



(12) 



by taking account of the (q— 1) dimensional vector struc- 
ture of the order parameter. The definition of Eq. ( jljj ) 
becomes a usual one for the scalar order parameter 
(q = 2). Since the prefactors of the L dependence in the 
finite-size scaling equations are canceled out, one may 
determine the critical temperature T c from the crossing 
point of the data of temperature dependence for differ- 
ent sizes as far as the corrections to finite-size scaling are 
negligible. We also plot the specific heat in Fig. 0(c). 
The specific heat peak becomes sharper when the system 
size is larger. From Figs. 0(a), 0(b) and 0(c), we find a 
clear phase transition at a finite temperature. 

The finite-size scaling plots of the order parameter, and 
the Binder parameter 

{M 2 )=L- 2p ' v f{[T-T c )L 1 ' v ^, (13) 
g = g({T-T c )L 1 l v ), (14) 

are given in Figs. 0(a) and 0(b), respectively. The esti- 
mated values for the critical temperature and the crit- 
ical exponents are T c = 1.222(4), l/z/= 1.52(4) and 
/3/i/=0.46(5). Here, the number in the parentheses de- 
notes the uncertainty in the last digit. The estimated 
values are consistent with the previous studies §-§, and 
the obtained exponents are close to those of the 3D XY 
model. We have given the estimate using the finite-size 
scaling analysis to check the consistency of the calcu- 
lation. Precisely speaking, due to small system sizes, 
our estimate of T c is a little bit lower than the accurate 
estimate |(],0], which results in a little bit smaller j3/u. 



To discuss more accurate estimates of T c and the criti- 
cal exponents, calculations with larger system sizes are 
preferable. 

It is not easy to estimate the specific heat exponent a 
from the specific heat data using the finite-size scaling 
analysis, if a is negative |6|. We may use the tempera- 
ture derivative of the specific heat, which is singular at T c 
even if a is negative. We plot the temperature derivative 
of the specific heat in Fig. 0. Since we directly calculate 
the density of states by using the Wang-Landau algo- 
rithm, it is easy to compute the temperature derivative 
of the specific heat from the moments of energy. Using 
the finite-size scaling relation 



M =L («-i)/» f ( {T _ Tc)L i/^ t 



(15) 



we estimate the exponent a as -0.04(6). Although our 
estimate has a relatively large error bar due to the small 
system size, our result suggests that a is negative. 

Next turn to the 4-state Potts model. We plot the 
temperature dependence of the order parameter and the 
Binder parameter in Fig. 0(a) and 0(b) , respectively. We 
also show the specific heat in Fig7^|(c) . We here make 
a comment on the statistical errors for the estimate of 
the density of states. They becomes larger for higher 
q and larger L. The errors for the specific heat curve 
for larger L are larger than the thickness of the curve 
in Fig. 0(c). From Fig. [| we find that there exits a 
clear finite-temperature phase transition. The finite-size 
scaling plots of the order parameter and the Binder pa- 
rameter are given in Figs. 0(a) and 0(b), respectively; 
we estimate T c , 1/u and (3/v as 0.669(4), 1.41(4) and 
0.44(6). They are compatible with the previous study 
ifjlfl , and the obtained exponents are close to those of 
the 3D Heisenberg model. Precisely, our estimate of T c 
is a little bit lower than the accurate estimate jl2); the 
situation is the same as the case of q = 3. 

The order parameter, the Binder parameter and the 
specific heat for the 5-state Potts model are given in 
Figs. (a), 0(b) and 0(c), respectively. There is no 
anomaly in the specific heat. There is no crossing in 
the Binder parameter at finite temperatures. The criti- 
cal temperature is down to zero. However, the value of 
the Binder parameter at T = is finite, and it may be- 
come constant for large enough L\ it is not clear whether 
T = is critical or not. To determine this point, more 
elaborate study with larger sizes are necessary. 

We show the order parameter, the Binder parame- 
ter and the specific heat for the 6-state Potts model in 
Figs. 0(a), 0(b) and 0(c), respectively. In this case, the 
maximum linear size is L — 14. There is no anomaly in 
the specific heat. The value of the Binder parameter at 
T = is very small. This means that the distribution 
of the order parameter is Gaussian; in other words, the 
system is disordered. We may conclude that there is no 
order even at T = 0. 
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IV. GROUND-STATE PROPERTIES 



S /N = 1.2717(1) + 0.17(4) x (l/N), q = 6, (19) 



In this section, we focus on the ground-state properties 
of the AF Potts models. First, we consider the ground- 
state entropy per spin, So/N; the entropy (in units of 
ks) is calculated through So = lng(0) 

In Fig. ||, we show the size dependence of the ground- 
state entropy per spin for the g-state AF Potts model on 
a simple cubic lattice. We plot Sq/N as a function of 
l/N, and find a linear l/N dependence. Using the least 
square method, for q = 3 we have 

S /N = 0.3670(1) + 1.97(4) x (l/N), q = 3, (16) 

where the number in the parentheses denotes the uncer- 
tainty in the last digit. Wang et al. Q estimated the 
ground-state entropy per spin (N — ► oo) as 0.3673 from 
the data of L— 4 and 8. Our estimate, 0.3670, is a lit- 
tle bit smaller than their estimate. The low-temperature 
phase of the AF 3-state Potts model is the BSS state. 
The typical configuration of the BSS state is that all the 
spins on the sublattice A take one of the three states, and 
the spins on the sublattice B take one of the other two 
states randomly. Then, the lower bound for the ground- 
state entropy becomes S /N > (ln2)/2 + In 6 xJl/N) = 
0.3466 + 1.79 x (l/N). Our estimate, Eq. @, is, of 
course, larger than the lower bound. That is, there are 
many configurations other than the typical configurations 
of the BSS state. 

From the linear 1/iV dependence for the AF 4-state 
Potts model, we have 

So/N = 0.7148(1) + 1.92(6) X (l/N), q = 4. (17) 

As far as we know, there has been no study on the es- 
timate of the ground-state entropy for q — 4 and higher 
q. For the low-temperature phase of the 4-state Potts 
model, the typical spin configuration is as follows: The 
spins on the sublattice A take two of four states ran- 
domly, whereas those on the sublattice B take the other 
two states. Then, the lower bound for the ground- 
state entropy becomes So/N > In 2 + In 6 x (l/N) — 
0.6931 + 1.79 x (l/N). Our estimate, Eq. (|l|), is again 
larger than the lower bound. 

From the size dependence of the ground-state entropy 
for the AF 5-state Potts model, which is also shown in 
Fig. ||, we have 

So/N = 0.9997(1) + 2.03(6) x (l/N), q = 5. (18) 

The main contribution to the low-temperature phase of 
the 5-state Potts model is as follows: The spins on the 
sublattice A take two of five states randomly, whereas 
those on the sublattice B take the other three states. 
Then, the lower bound for the ground-state entropy be- 
comes S /N > (In 6)/2 + 2 In 5 x (l/N) = 0.8959 + 3.00 x 
(l/N). Our estimate, Eq. (|l8|), is again larger than the 
lower bound. 

Finally, for q = 6 we have 



by using the least square method. The main contri- 
bution of the 6-state Potts model is as follows: The 
spins on the sublattice A take three of six states ran- 
domly, whereas those on the sublattice B take the other 
three states. Our estimate, Eq. (|l9|), is again larger than 
the lower bound for the ground-state entropy; So/N > 
In3 + 2hi5 x (l/N) = 1.0986 + 3.00 x (1/iV). 

It is interesting to note that the ground-state entropy 
has a similar l/N dependence for all q — 3,4, 5, although 
there occurs the finite-temperature phase transition for 
<7 = 3,4 and the zero-temperature phase transition for 
g = 5. The size dependence of q = 6 is similar to that for 
q = 3, 4, 5. But the slope of l/N dependence for q = 6 is 
smaller than that of others. It is not clear whether this 
small difference is related to the existence of the antifcr- 
romagnetic order. 

We find that the ground-state entropy is larger than 
the contribution from the typical configurations of the 
BSS state. In order to look into the ground-state proper- 
ties more carefully, let us consider the proportion of vis- 
iting the typical BSS ground states among all the ground 
states. As an example, we treat the case of q = 3. Since 
this proportion becomes very small for larger sizes, we 
have checked it for a smaller system, that is, L = 6. The 
estimate of the ground state entropy per spin So/N for 
L = 6 is 0.376. Thus, the expected value for the pro- 
portion of the typical BSS ground states among all the 
ground states is 6 • 2W 2 ) / e °- 376Ar = 0.010 (AT=216). Ac- 
tually, the typical BSS ground states were visited 0.010 
times as frequently as all the ground states, which is con- 
sistent with the theoretical expectation. In other words, 
we sample ground states uniformly. 

Next we study the distribution of the ground states. 
As an example, we again deal with the case of q = 3. 
Let us denote the number of each component per spin as 
nf' B (z = 1,2,3); then, 

0<n,<l/2 and rrH=nf-nf, (20) 

where the staggered magnetization vrii is given in Eq. (0). 

A B — 

We show the distribution function of n i ' for L—8, 10, 
12, 14 and 16 in Fig. ^. The same types of lines for differ- 
ent sizes are used as in Fig. |l|. The distribution is sharper 
for larger sizes. In plotting the data in Fig. ||, we have 
chosen the component 1 and the sublattice A such that 
is the largest. The distribution function P(nf' B ) is 
normalized such that 

/ P(4 > B ) dnt B 

is independent of the size. For the typical BSS states, 
we expect the (5-function distribution for the sublattice 
A, P(n?) = 5(nf - 1/2), P(n^) = 6(n£ 3 ); P(n B ) + 
P(n B ) becomes the Gaussian distribution around 1/4 
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and P(nf) — S(nf) for the sublattice B. Our results 
shown in Fig. ^| are close to those of the typical BSS 
states, but we can see a clear deviation from the typi- 
cal BSS states and the size dependence. In other words, 
there are many configurations other than the typical con- 
figurations of the BSS state. Similar behavior is also ob- 
tained for the ground states of higher-state Potts models. 
The deviation from the typical BSS states becomes larger 
for higher q. 



V. SUMMARY AND DISCUSSIONS 



To summarize, we have applied a newly proposed 
Monte Carlo algorithm, the Wang-Landau algorithm 
H|, to the study of the 3D AF g-state Potts models. 
We obtain the finite-temperature phase transition for q= 
3 and 4, whereas the transition temperature is down to 
zero for q=5. For q—6 there exists no order for all the 
temperatures. We also study the ground-state properties. 
From the analysis of the size-dependence of the ground- 
state entropy, we find that the ground-state entropy is 
larger than the contribution from the typical configura- 
tions of the BSS state for q = 3. The same situations are 
found for q = 4, 5 and 6. 

We have confirmed again the efficiency of the Wang- 
Landau algorithm. For the study of only the critical phe- 
nomena near the critical point, other methods may have 
advantage. However, in order to make the systematic 
study for all the energy space, we can treat larger sys- 
tems by using the Wang-Landau algorithm, especially for 
higher g-state Potts models. 

In the present paper we have studied the ground-state 
entropy for the 3D AF Potts models. We had better men- 
tion that the accurate estimate of the nonzero ground- 
state entropy is easily obtained for the 2D AF Potts mod- 
els, of course; the obtained data are consistent with the 
previous studies p3[ . It is also interesting to apply the 
Wang-Landau algorithm to the systematic study of both 
the ground-state properties and the phase transitions for 
more complicated systems, such as spin glass problems. 
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(T - T c ) L 

FIG. 2. Scaling plots of the order parameter (a) and the 
Binder parameter (b) for the 3D AF 3-state Potts model. 
The linear system sizes are L = 8, 10, 12, 14 and 16. The 
estimated values for the critical temperature and the critical 
exponents are T c = 1.222, l/v= 1.52 and /3/V=0.46; 
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FIG. 1. Temperature dependence of the squared order pa- 
rameter (a), the Binder parameter (b), and the specific heat 
(c) for the 3D AF 3-state Potts model. The linear system 
sizes are L — 8, 10, 12, 14 and 16. The dot-dashed, dashed, 
short-dashed, dotted and solid lines are used for L=8, 10, 12, 
14 and 16, respectively. The temperature is represented in 
units of J/ks- 



Oil— 



-5 



-10 





I I I 1 1 


i i 

q=3; 




.^• •^ 7 X \'i 






|\ 
\ 






L=16\ 


'- 











1.25 



1.5 



FIG. 3. Temperature derivative of the specific heat for the 
3D AF 3-state Potts model. The linear system sizes are L = 
8, 10, 12, 14 and 16. 
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FIG. 5. Scaling plots of the order parameter (a) and the 
Binder parameter (b) for the 3D AF 4-state Potts model. 
The linear system sizes are L = 8, 10, 12, 14 and 16. The 
estimated values for the critical temperature and the critical 
exponents are T c = 0.669, l/v= 1.41 and /3/u=0A4; 



FIG. 4. Temperature dependence of the squared order pa- 
rameter (a), the Binder parameter (b), and the specific heat 
(c) for the 3D AF 4-state Potts model. The linear system sizes 
are L = 8, 10, 12, 14 and 16. The temperature is represented 
in units of J/ks- 
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FIG. 6. Temperature dependence of the squared order pa- 
rameter (a), the Binder parameter (b), and the specific heat 
(c) for the 3D AF 5-state Potts model. The linear system sizes 
are L = 8, 10, 12, 14 and 16. The temperature is represented 
in units of J/fcs. 
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FIG. 7. Temperature dependence of the squared order pa- 
rameter (a), the Binder parameter (b), and the specific heat 
(c) for the 3D AF 6-state Potts model. The linear system sizes 
are L = 8, 10, 12 and 14. The temperature is represented in 
units of J/fcs. 
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FIG. 8. Size dependence of the ground-state entropy per 
spin for the 3D AF g-state Potts models; g=3, 4, 5, and 6. 
The linear system sizes are L = 8, 10, 12, 14 and 16 for g=3 
to 5 and L = 8, 10, 12 and 14 for g=6; N = L 3 . 
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FIG. 9. Distribution function P(x) for the ground states 
of the 3D AF 3-state Potts models. Here x stands for 
n t' B (* ~ 1, 2 9 3), the number of each component per spin 
on the sublattice A or B. The linear system sizes are L = 8, 
10, 12, 14 and 16. The same types of lines for different sizes 
are used as in Fig. |l[ 
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